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We present a first-principles study of the atomic, electronic and magnetic properties of two dimen- 
sional (2D), single and bilayer ZnO in honeycomb structure and its armchair and zigzag nanoribbons. 
In order to reveal the dimensionality effects, our study includes also bulk ZnO in wurtzite, zincblende 
and hexagonal structures. The stability of 2D ZnO, its nanoribbons and flakes are analyzed by 
phonon frequency, as well as by finite temperature ab-initio molecular dynamics calculations. 2D 
ZnO in honeycomb structure and its armchair nanoribbons are nonmagnetic semiconductors, but 
acquire net magnetic moment upon the creation of zinc vacancy defect. Zigzag ZnO nanoribbons are 
ferromagnetic metals with spins localized at the oxygen atoms at the edges and have high spin po- 
larization at the Fermi level. However, they change to nonmagnetic metal upon termination of their 
edges with hydrogen atoms. From the phonon calculations, the fourth acoustical mode specified as 
twisting mode is also revealed for armchair nanoribbon. Under tensile stress the nanoribbons are 
deformed elastically maintaining honeycomb like structure, but yield at high strains. Beyond yield- 
ing point honeycomb like structure undergo a structural change and deform plastically by forming 
large polygons. The variation of the electronic and magnetic properties of these nanoribbons have 
been examined under strain. It appears that plastically deformed nanoribbons may offer a new class 
of materials with diverse properties. 

PACS numbers: 73.22.-f, 75.75.+a, 63.22.-m 



I. INTRODUCTION 

Graphene, a monolayer layer of carbon atoms in hon- 
eycomb structure, is offering exceptional properties"^ 
which may lead to important applications in various 
fields. Normally, two dimensional (2D) graphene is 
scmimetallic and its electrons and holes behave like 
a massless Dirac fermion. Whereas 2D boron-nitride 
(BN)J3 Group III-V analogue of graphene, in ionic hon- 
eycomb structure is a wide band gap semiconductor. Un- 
usual properties of graphene and BN nanoribbons have 
been revealed extensively in recent papersP"^ More re- 
cently, based on state-of-the art first principles calcula- 
tions it was predicted that Si and Ge^l, even binary com- 
pounds of Group IV elements and III-V compounds^ 
can form 2D stable monolayer honeycomb structures. 
Earlier studies on ZnO and its nanowir es gav e first in- 
dications that graphitic ZnO can existsfSED Very thin 
nanosheetsp^ nanobeltsp^ nanotubes^ and nanowires^ 
of ZnO have already been synthesized. Two monolayer 
thick ZnO(OOOl) films grown on Ag(lll) were reported.^ 

Because of its wide band gap of ^3.3 eV and large 
exciton binding energy of 60 meV leading to vast op- 
toelectronic applications including light-emitting diodes 
and solar cells 123221 ZnO has been the subject of various 
researches. It was reported that the vacancy defects can 
be intentionally created by electron irradiation method.^ 
It has been also reported that Zn vacancy induces fer- 
romagnetism in ZnO thin films and nanowires without 
any need of doping with transition-metal atoms These 
magnetic properties might provide superior advantages 
in biomedical applications because of non-toxic nature of 
ZnO as opposed to transition metal ions. 



In this paper a comprehensive study of the atomic, 
electronic and magnetic properties of monolayer, bilayer 
and nanoribbons of II- VI ionic ZnO compound in hon- 
eycomb structures are carried out using first-principles 
calculations. In order to reveal the dimensionality ef- 
fects, we started with the energetics and electronic energy 
bands of ZnO in different 3D bulk crystalline structures 
and compared them with those of 2D ZnO honeycomb 
structure. Our analysis based on phonon dispersions and 
finite temperature ab-initio molecular dynamics calcula- 
tions provides evidence for the stability of free-standing 
2D monolayer, bilayer and quasi ID nanoribbons of ZnO 
honeycomb structures. These structures can be in dif- 
ferent local minima on the Born-Oppenheimer surface, 
in spite of the fact that they are not synthesized yet. 
We found that 2D monolayer and bilayer ZnO are non- 
magnetic semiconductors, but attain magnetic proper- 
ties upon creation of Zn- vacancy defect. ZnO nanorib- 
bons exhibit interesting electronic and magnetic proper- 
ties depending on their orientation. While armchair ZnO 
nanoribbons are nonmagnetic semiconductors with band 
gaps varying with their widths, bare zigzag nanoribbons 
are ferromagnetic metals. These electronic and mag- 
netic properties show dramatic changes under elastic and 
plastic deformation. Hence, ZnO nanoribbons can be 
functionalized by plastic deformation. Results obtained 
in this study indicates that 2D monolayer and bilayer 
ZnO honeycomb structures and quasi ID armchair and 
zigzag nanoribbons display unusual electronic, magnetic 
and mechanical properties, which hold the promise of in- 
teresting technological applications. 



2 



II. MODEL AND METHODOLOGY 

We have performed first-principles plane wave calcu- 
lations within Density Functional Theory (DFT) using 
PAW potentials.^ The exchange correlation potential 
has been approximated by Generalized Gradient Approx- 
imation (GGA) using PW9lP functional both for spm- 
polarized and spin-unpolarized cases. Recently, spin- 
polarized calculations within DFT have been used suc- 
cessfully to investigate magnetic properties of vacancy 
defects in 2D honeycomb structures. Also interesting 
spintronic properties of nanoribbons have been revealed 
using spin-polarized DFT P The success of spin-polarized 
DFT calculations has been discussed in Ref2"Tl 

All structures have been treated within supercell ge- 
ometry using the periodic boundary conditions. A plane- 
wave basis set with kinetic energy cutoff of 500 eV has 
been used. The interaction between ZnO monolayers in 
adjacent supercells is examined as a function of their 
spacing. Since the total energy per cell has changed less 
than 1 meV upon increasing the spacing from 10 A to 
15 A , we used the spacing of ~ 10 A in the calcula- 
tions. In the self-consistent potential and total energy 
calculations the Brillouin zone (BZ) is sampled by, re- 
spectively (15x15x15), (25x25x1) and (25x1x1) special k- 
points for 3D bulk, 2D honeycomb and ID (nanoribbons) 
ZnO. This sampling is scaled according to the size of su- 
perlattices. For example, BZ was sampled by (3x3x1) 
special k-points for defect calculations using (7x7) su- 
percell of 2D ZnO honeycomb structure. All atomic po- 
sitions and lattice constants are optimized by using the 
conjugate gradient method, where the total energy and 
atomic forces are minimized. The convergence for energy 
is chosen as 10 -5 eV between two steps, and the maxi- 
mum Hellmann-Feynman forces acting on each atom is 
less than 0.02 eV/A upon ionic relaxation. The pseu- 
dopotentials having 12 and 6 valence electrons for the Zn 
(Zn: 4s 2 3d 10 ) and O ions (O: 2s 2 2p 4 ) are used. Numeri- 
cal plan e wav e calculations are performed by using VASP 
package! 28 * 29 ^ While all numerical calculations of struc- 
ture optimization, electronic energy and pho non di sper- 
sions are carried out within GGA using VASPp^EH some 
of the calculations are checked also by using PWSCF 30 
and SIESTA 31 softwares. Therefore all pertaining discus- 
sions are based on the results obtained by using VASP 
software unless it is stated otherwise. 

Since DFT within GGA underestimates the band gap, 
frequency-dependent GWq calculations are carried outP^l 
Screened Coulomb potential, W, is kept fixed to initial 
DFT value Wo and Green's function, G, is iterated five 
times. Various tests regarding vacuum separation, ki- 
netic energy cut-off potential, number of bands, k-points 
and grid points are made. Final results of GWq cor- 
rections are obtained using (12x12x1) k-points in BZ, 
20 A vacuum separation, 400 eV cut-off potential, 160 
bands and 64 grid points. In addition to GWq, we per- 
formed also GW calculations in order to make compari- 
son with earlier available studies. While GWq corrections 



are successfully applied to 3D and 2D ZnO, its applica- 
tion to quasi ID nanoribbons is hindered by large number 
of atoms. 



III. 3D BULK AND 2D HONEYCOMB ZNO 
CRYSTAL 

We first consider 3D bulk ZnO, which are in wurtzite 
(wz-ZnO), zincblende (zb-ZnO) and hexagonal (h-ZnO, 
or graphite-like) crystals. Atoms in wz- and zb-ZnO 
structures are four fold coordinated through tetrahe- 
drally directed suborbitals, whereas the atoms in h-ZnO 
crystal are three fold coordinated through sp 2 -orbitals. 
Wurtzite structure is found to be the thcrmodynamically 
most stable phase of ZnO. The cohesive energy per Zn-0 
pair is calculated by using the expression, 

E c = E T [ZnO] -E T [Zn] - E T [0] (1) 

in terms of the total energy of the optimized crystal 
structure of ZnO, Et[ZuO] per Zn-0 pair, the total en- 
ergies of free Zn and O atoms Et\Zv\ and Et[0]. The 
equilibrium cohesive energies of wz-ZnO, zb-ZnO and h- 
ZnO structures are found to be 8.934, 8.919 and 8.802 
eV per Zn-0 pair, respectively. For wz-ZnO crystal, the 
hexagonal lattice constants of the optimized structure in 
equilibrium are a = 3.280 A, c/a = 1.616. The devia- 
tion of c/a from the ideal value of 1.633 imposes a slight 
anisotropy in the lengths of tetrahedrally directed Zn-0 
bonds. While the length of three short bonds is 2.001 A, 
the fourth bond is slightly longer and has the length of 
2.007 A. The zincblende structure in Td symmetry has 
cubic lattice constant a — 3.266 A and four tetrahe- 
drally coordinated Zn-0 bonds having uniform length, 
d = 2.001 A . The h-ZnO structure has hexagonal lattice 
constants a = 3.448 A, c/a = 1.336 and d = 1.990. The 
structural parameters of these three bulk ZnO crystals 
are shown in Fig. [T] The lattice constant of wz-ZnO, a 
(c/a) is measured between 3.247 (1.6035) A and 3.250 
(1.602) A using different methods.^ The calculations 
based on ab-initio LCAO method with all-electron Gaus- 
sian type basis set predict a=3.286 A and c/a=1.595.^ 
All our results related with the structural parameters are 
in good agreement with the experimental and t heoretical 
values within the average error of ~ 1 %! 33 l 34 l 

The optimized atomic structure and corresponding 
electronic band structure of 3D ZnO crystals are pre- 
sented in Fig. [T] All wz-ZnO, zb-ZnO and h-ZnO crys- 
tals are direct band gap semiconductors with calculated 
band gaps being Ec=0-75, 0.65 and 0.96 eV, respectively. 
Highest valence band has 0-2p character; the states of 
the lowest conduction band is formed from Zn-^p and 
Zn-^s orbitals. Valence band consists of two parts sepa- 
rated by a wide intra band gap. The lower part at ~ -18 
eV is projected mainly to 0-2s orbitals. The upper part 
is due to mainly Zn-3ci and 0-2s orbitals. The differ- 
ences in the band structure of different three 3D crystals 
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Figure 1: (Color online) Cohesive energies Ec per Zn-0 pair, 
band gaps Eg, Zn-O bonds, atomic and energy band struc- 
tures of 3D wurtzite (a), zincblende (b), and hexagonal (c) 
crystals of bulk ZnO. Red-small and blue-large balls corre- 
spond to O and Zn atoms, respectively. The gap between 
valence and conduction bands are shaded and the zero of en- 
ergy is set at the Fermi level Ef- All structures including 
lattice constants are fully optimized. 




Figure 2: (Color online) Phonon dispersion curves calculated 
by force constant method for 2D monolayer ZnO. Atomic 
structure and hexagonal lattice vectors are shown by inset. 



zb-ZnOp3 AH electron LAPW calculations predicted the 
band gap of wz-ZnO 0.77 eV using LDA. which is cor- 
rected by GW calculations to 2.44 eV.^I 

Charge transfer from Zn atoms to O atoms is a mea- 
sure of the ionicity of ZnO crystal. We calculated the 
amount of charge on constituent Zn and O atoms in 3D 
crystals by performing the Lowdin analysis^ in terms 
of the projection of plane-waves into atomic orbitals. By 
subtracting the valencies of free Zn and O atoms from the 
calculated charge values on the same atoms in 3D crys- 
tals we obtain the charge transfer Sq from Zn to O. The 
calculated value of charge transfer for both wz-ZnO, zb- 
ZnO and h-ZnO is found to be 1.41 electrons. The values 
of charge transfer calculated with the Bader analysis^! 
are Sq — 1.22, Sq = 1.17 and Sq = 1.20 electrons for 
wz-ZnO, zb-ZnO and h-ZnO, respectively. Same analy- 
sis performed with local basis set using SIESTA^H yields 
significantly lower values of charge transfer, 5q = 0.90, 
8q = 0.88 and Sq = 0.91 electrons for wz-ZnO, zb-ZnO 
and h-ZnO, respectively. This analysis clearly indicates 
that a significant amount of charge is transferred from 
low electronegative Zn atom to high electronegative O 
atom. However the values of Sq may scatter owing to 
the ambiguities in placing boundary between Zn and O 
in crystalline structure. 



IV. 2D ZNO HONEYCOMB STRUCTURE 



become pronounced in the lower part of the conduction 
band. It should be noted that the band gaps of bulk ZnO 
is highly underestimated by DFT calculations.^ The ex- 
perimentally measured band gap of wz-ZnO is ^3.37 eV 
at room temperature.^ We performed GWq calculations 
to correct the band gaps calculated within GGA. Our 
results for wz-ZnO, zb-ZnO and h-ZnO are, respectively, 
3.29, 3.04 and 3.32 eV. As for earlier studies, calculations 
with GW corrections reported a band gap of 3.59 eV for 



The structure of monolayer ZnO is optimized using pe- 
riodically repeating supercell having hexagonal lattice in 
2D and the spacing of 10 A between ZnO planes. The op- 
timized structure was planar and the magnitudes of the 
Bravais vectors of the hexagonal lattice are found to be 
a x =a 2 = 3.283 A, and the Zn-0 bond, d =1.895 A (see 
Fig. [2]). In contrast to 2D puckered honeycomb structures 
of Si, Ge and compounds of III and V group elements 
both lying below the first row, any honeycomb structure 
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Figure 3: (Color online) Electronic structure of the 2D mono- 
layer of ZnO. (a) Brillouin zone corresponding to 2D hexago- 
nal lattice and difference charge density, Ap. (b) Energy band 
structure and density of states projected to the selected or- 
bitals (PDOS) of constituent atoms. The band gap is shaded 
(painted in yellow) and the zero of energy is set at the Fermi 
level. Band edges corrected by GWo are indicated by small ( 
light blue) circles. The band gap enlarged after GWo correc- 
tion is shaded (painted light blue). 



including one element from the first row form planar hon- 
eycomb structure like graphene, BN and SiCP^2D mono- 
layer of ZnO forming planar honeycomb structure is not 
an exception. The calculated structural parameters are 
significantly larger than those of graphene and BN hon- 
eycomb structure due to fact that Zn has larger radius^ 
than that of B,C,N and O atoms. The planar structure 
of 2D ZnO is tested by displacing Zn and O atoms ar- 
bitrarily from their equilibrium positions by 0.5 A and 
subsequently by re-optimizing the structure. Upon opti- 
mization the displaced atoms have recovered their origi- 
nal positions in the same plane. It should be noted that, 
the length of Zn-0 bonds of 2D ZnO honeycomb struc- 
ture is smaller than that in the 3D bulk (wz, zb) crystals, 
since sp 2 bonding in the former is stronger than the tetra- 
hedrally coordinated sp 3 bonding in the latter. Similar 
trend is also found in C, BN and SiC honeycomb struc- 
tures. The interaction between ZnO planes appears to 
slightly weaken the Zn-0 bonds of h-ZnO crystal. As a 
result the length of the Zn-0 bonds becomes larger than 
that of 2D ZnO honeycomb structure. The cohesive en- 
ergy of 2D monolayer of ZnO is calculated to be 8.419 eV 
per Zn-0 pair. This energy is 0.5 eV smaller than that 
of 3D wz-ZnO. 



A. Atomic structure and stability analysis 

It should be noted that 2D honeycomb structure deter- 
mined by the structure optimization using conjugate gra- 
dient method may not be necessarily stable. One needs 
more stringent tests to assure the stability of this truly 
2D structure. As a matter of fact, it was argued earlier 
that any crystal structures con sisting of a truly 2D mono- 
layer crystals cannot be stableP^U Here, we summarize 
our analysis on the stability of planar 2D hexagonal struc- 
ture of ZnO using calculation of phonon dispersion curves 
as well as ab-initio finite temperature molecular dynam- 
ics calculations. 

We calculated optical and acoustical branches of 
phonon frequency using the density functional pertur- 
bation theory (DFPT) with plane wave method as im- 
plemented in PWSCF software and the force constant 
method^ with forces calculated using VASP. The former 
method includes the polarization effects so that the split- 
ting of longitudinal and transverse optical modes at the 
T-point (namely LO-TO splitting) is obtained. The force 
constant method cannot yield the LO-TO splitting and 
is also very sensitive to the mesh size in calculating forces 
under atomic displacements and hence in setting up the 
dynamical matrix. In fact, the imaginary frequencies of 
ZA branch (so called out of plane acoustical branch) oc- 
cur around the T-point as an artifact of insufficient mesh 
size. However, all imaginary frequencies around the T- 
point (corresponding to an instability to be induced by 
acoustical waves with large A) are removed by using very 
fine mesh. Earlier, imaginary frequencies of ZA modes 
near the T-point calculated for 2D honeycomb structures 
of Ge and III-V compounds inducing similar fortuitous 
instabilities for large A were also removed by using finer 
meshP^ 

In Fig. [2] we present phonon dispersion calculated for 
2D monolayer of ZnO honeycomb structure using force 
constant method, where atomic forces are calculated by 
VASP. Calculated phonon modes all being positive across 
the BZ strongly support that 2D monolayer of ZnO is a 
stable structure corresponding to a local minimum on 
the Born-Oppcnheimer surface. Our phonon dispersion 
curves are in agreement with those calculated in Ref|44| . 

Furthermore, we have tested the stability of 2D ZnO 
monolayer using finite temperature ab initio molecular 
dynamics (MD) calculations with time steps of 2xl0~ 15 
seconds. To lift the constraints to be imposed by a 
small unit cell, we considered a (7x7) supercell of the 
2D hexagonal ZnO and raised its temperature from K 
to 750 K in 0.1 ps. Then, we have kept the temperature 
of the system around 750 K for 2.5 ps. While all these 
calculations resulted in minor deformations, the honey- 
comb structure was not destroyed. It should be noted 
that these calculations may not be conclusive, since 2.5 ps 
cannot be sufficient to represent all the statistics, but this 
picture is the best one can see from the exiguous window 
limited by the computational cost imposed by ab-initio 
MD method. In the following sections we will present 
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additional arguments related with electronic structure, 
which further corroborate the stability of 2D ZnO struc- 
tures. 



B. Electronic Structure 

The difference charge density and the electronic en- 
ergy bands together with the partial density of states are 
presented in Fig. [3] Contour plots of total charge indi- 
cate high density around O atoms. The difference charge 
density is calculated by subtracting charge densities of 
free Zn and O atoms from the total charge density of 2D 
ZnO, i.e. Ap = p Zn o - Pzn - Po- High density con- 
tour plots around O atoms protruding towards the Zn-0 
bonds indicate charge transfer from Zn to O atoms. This 
way the Zn-0 bond acquires an ionic character. The 
charge transfer from Zn to O, 5q is analyzed by using 
different schemes. The charge transfer values calculated 
by Lowdin method using PWSCF , Badei 39 analysis us- 
ing VASP, and local basis set analysis using SIESTA are, 
respectively, 5q=l.35, 1.18 and 0.87 electrons. Interest- 
ingly, as compared to 3D bulk h-ZnO crystal, the charge 
transfer from cation to anion of 2D monolayer ZnO is 
slightly decreased. This is due to the change from sp 3 
hybrid orbital in h-ZnO to sp 2 hybrid orbital in honey- 
comb structure. 

Two dimensional ZnO is a direct band gap semicon- 
ductor with a gap value of f .68 eV. However, the actual 
band gap is expected to be larger. The bands are cor- 
rected using GWq method^ and the direct band gap at 
the T-point of BZ is found to be 5.64 eV. Much recently, 
it has been reported that the band gap of 2D- ZnO is 
calculated to be 3.57 eV with GW corrections.^ Sim- 
ilar to the bands of 3D bulk ZnO crystals, the upper 
part of the upper valance band are derived mainly from 
0-2p orbitals, whereas the lower part has Zn-3d charac- 
ter. The bands at the edges of conduction and valence 
bands along the T-K direction are derived from bond- 
ing and antibonding combination of 0-2p z and Zn-4p 2 
orbitals forming tt- and 7r* -states. The highest valence 
band along the T-K direction has mainly 0-2p z but small 
Zn-4p z orbital contribution, while the lowest conduction 
band is composed mainly from Zn-4p z but small 0-2p z 
orbitals. Small Zn-4p z contribution is also confirmed by 
PDOS. Consequently, planar geometry of 2D monolayer 
of ZnO is expected to be attained by the bonding com- 
bination of p z -states. The same situation occurs for 2D 
monolayer of BN honeycomb structure, which is an ionic 
lll-V compound with wide band gap between tt- and tt*- 
states. The planar stability graphene is also maintained 
by 7r-bonds. 



C. Vacancy Defects and Antisite 

It was shown that vacancies have remarkable effects on 
electronic and magnetic properties of 2D graphene and 
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Figure 4: (Color online) Vacancy defects in a (7x7) supercell 
of the monolayer of ZnO. (a) Relaxed atomic structure around 
the Zn-vacancy with isosurfaces of the difference charge den- 
sity of spin-up and spin-down states, Ap = p(f) = p(i). En- 
ergies Ei of localized states measured from the top of the 
valance band and the net magnetic moment per supercell in 
Bohr Magneton ps- Blue-large and red-small balls indicate 
Zn and O atoms indicated by arrows, (b) O-vacancy. (c) 
Zn+O divacancy. (d) O+Zn antisite where neighboring O 
and Zn atoms are exchanged. The (7x7) supercells each in- 
cluding one type of various vacancy defects are delineated by 
solid lines. 



graphene nanoribbonsP^HSID To the best of our knowl- 
edge, the effects of vacancies in monolayer of ZnO have 
not been treated yet. It was reported that the vacancy 
defects can be intentionally created by electron irradi- 
ation method^ on ZnO thin films. We investigated the 
effects of Zn-, 0-, Zn+O-di vacancy and O+Zn-antisite in 
a periodically repeating (5x5) as well as (7x7) supercells. 
The vacancy-vacancy coupling in the larger (7x7) super- 
cell was reduced significantly. Flat bands in the band 
gap have charge density localized at the defect site. The 
width of these flat bands can be taken as the measure 
of the strength of the direct and indirect (via the hop- 
ping through the orbitals in the lattice) vacancy-vacancy 
coupling. The largest width of such a band is already 
small (less than 50 meV). The average energy of these flat 
bands from the top of the valence band can be taken as 
the localized states of an individual vacancy defect. Our 
results obtained using a (7x7) supercell are presented in 
Fig. [4] and discussed in the rest of the section. 

The GGA band gap of a defect-free, 2D monolayer ZnO 
increases by ^0.08 eV in a (7x7) supercell including a 
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vacancy defect. This is an artifact of the model, which 
mimics an individual vacancy in supercell having limited 
size. First we generated a Zn-vacancy by removing a sin- 
gle Zn atom from the monolayer of ZnO in honeycomb 
structure as shown in Fig.[4](a). Upon the structure opti- 
mization, two O atoms around vacancy are departed from 
the plane. Similar to the vacancies in graphene and BN, 
Zn-vacancy induces a local magnetization in the system. 
Isovalue surfaces of spin density difference Ap clearly re- 
veals the magnetism around vacancy. The calculated to- 
tal magnetic moment is 2 pg per unit cell. Note that 
similar Zn-vacancy in the 3D bulk ZnO crystal induces 
a magnetic moment of ~ 1 ps- It should be noted that 
the magnetic moment of Zn-vacancy arises from the un- 
paired electrons of oxygen atoms around the Zn vacancy. 
Therefore it is important to determine precisely whether 
the vacancy induced localized states are occupied. The 
positions of these localized states of 2D ZnO calculated 
for the (7x7) supercell and their occupancy unambigu- 
ously comply with the calculated magnetic moment . The 
Zn-vacancy in a repeating (7x7) supercell also modifies 
the electronic structure. The value of the band gap of 
defect-free ZnO changes from 1.68 eV to 1.75 eV, and 
a spin up localized state appear 0.25 eV above the top 
of the valence band. Figure [4] (b) presents our results 
for O-vacancy. Unlike the case of Zn-vacancy, the mono- 
layer of ZnO containing an O-vacancy is nonmagnetic. 
Zn atoms around vacancy with coordination number 2 
prefer to stay in the same plane of the other atoms and do 
not induce any magnetization. The band gap is slightly 
modified to 1.80 eV. As for Zn+O divacancy in Fig.|4](c), 
it is again nonmagnetic. The band gap is also modified 
to 1.46 eV. Four occupied states associated with diva- 
cancy occur as localized states in the band gap. Finally, 
we consider the antisite defect. The resulting relaxed 
structure is given in Fig.|4](d). The antisite O is pushed 
away from Zn atoms and makes bonds with a nearby 
O atom. The system does not show any magnetization. 
The energy of final structure with antisite defect is ~ 5 
eV higher (i.e. energetically less favorable) that of the 
defect-free 2D monolayer of ZnO. The antisite induces 
two localized states in the band indicated in the figure. 
Also it is noted that the magnetic moments calculated 
for single Zn- and O-vacancy do not agree with Lieb's 
theorem,^ which predicts the amount of magnetic mo- 
ments for carbon vacancies in 2D graphene. According 
to Lieb's theorem, the net magnetic moment per cell is 
expected to be p=l pb for Zn and O vacancies in Fig. [4] 
This might be related with the ionic bonding which is 
different from graphene or the existence of rf-orbitals in 
Zn. On the other hand, zero magnetization for Zn+O 
divacancy in Fig. [I] (c) is in agreement with the theorem. 



D. 2D ZnO bilayer 

Recently, Tusche et alW^ revealed two monolayer of 
ZnO grown on Ag(lll) substrate using surface x-ray 



diffraction and scanning tunneling microscopy. They also 
showed that the transition to the bulk wz-ZnO structure 
occurs in the 3-4 monolayer coverage. The Zn-0 bond 
length of the planar hexagonal structure measured 1.93 
A is slightly larger than the the value of 1.895 A calcu- 
lated for the the bond length of the monolayer of ZnO. 
This situation implies that the effect of the Ag(lll) sub- 
strate may be negligiblePS 

We investigated the atomic structure and stability of 
bilayer ZnO honeycomb structure. To determine the min- 
imum energy configuration we used periodic supercell ge- 
ometry and considered five different stacking types which 
are T±, T 2 , Hi, H2 and H3. In T\ configuration, Zn (O) 
atoms of the second layer are on top of the Zn (O) atoms 
of the first layer. T 2 is similar to Xi except that Zn atoms 
are above O atoms. Hi, H2 and H3 configurations are 
obtained by shifting one of the layers of Ti and T 2 on the 
other layer, so that Zn or O atoms of the second layer 
are placed above the center of the hexagons in the first 
layer. It is Hi (H2) if O (Zn) atoms of both layers face 
the centers of hexagons. H 3 corresponds to the configu- 
ration where O (Zn) atom of the first (second) layer face 
the center of second (first) layer. 

Among these five configurations we found that T 2 is 
energetically most favorable. H3, H2, Hi and T\ configu- 
rations have respectively 0.213, 0.312, 0.320 and 0.321 eV 
(per primitive cell) higher energies than T2 configuration. 
The length of Zn-0 bonds as well as the charge transfer 
from Zn to O in T2 configuration is slightly larger than 
the value of 1.895 A calculated for the the bond length 
of the monolayer of ZnO. This trend complies with the 
above discussion that the bond length increases by going 
from 2D monolayer to 3D bulk and implies that the effect 
of the Ag(lll) substrate may be negligible.^ 

The bilayer binding energy among two ZnO layers is 
calculated to be 0.162 eV (per Zn-0 pair) by subtracting 
the total energies of two individual monolayers from the 
total energy of bilayer ZnO in T2 configuration. Fur- 
thermore the layer-layer separations are calculated as 
4.02, 2.36, 3.80, 3.78 and 2.68 A for T u T 2 , H u H 2 and 
i?3 respectively. Hence, bilayer formation is exothermic. 
The calculations with LDA, which accommodate van der 
Waals interaction between layers better, yields 0.297 eV 
(per Zn-0 pair) binding energy between two ZnO lay- 
ers and 2.267 A layer-layer separation in T2 structure. 
The layer-layer separation for T 2 is experimentally re- 
ported as 2.314 AP3 The difference between the lattice 
constants of the monolayer 2D ZnO honeycomb structure 
and the ZnO bilayer is small. Owing to the relatively 
weak coupling between two ZnO monolayers, the calcu- 
lated electronic structure is similar to that of single layer 
ZnO honeycomb structure, except that the band gap de- 
creases to 1.44 eV for T 2 . This band gap increases to 5.10 
eV after the GWq correction, which is still smaller than 
the GWo corrected band gap of 2D monolayer ZnO. 

We next address the question whether the structure 
of 2D bilayer formed by the T 2 stacking of the ZnO bi- 
layer is stable is examined by the finite temperature ab- 
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Figure 5: (Color online) (a) Phonon dispersion curves of 
a-ZnO(9)calculated by force constant method. Right panel 
presents the magnified view of the low lying modes. Here the 
well known out of plane, transverse and longitudinal acousti- 
cal modes are labeled, respectively, as ZA, TA and LA. The 
fourth acoustical twisting mode, which is indigenous to quasi 
ID nanoribbons, is labeled as TW, because this mode cor- 
responds to the twisting of the nanoribbon along the axis 
passing through the middle of the ribbon in the infinite di- 
rection, (b) Density of states (DOS) calculated for 2D planar 
ZnO (shown by solid lines); armchair a-ZnO(9) (shown by 
dashed lines); a finite size flake, i.e. a 5 unit cell long a- 
ZnO(9)(shown by dotted lines). The discrete frequencies due 
to finite size structure are gaussian broaden with <r=8 cm -1 in 
order to compare with continuous DOS of 2D layer and quasi 
ID ribbon. Quasi ID (nanoribbon) and OD (flake) honeycomb 
structures are shown by inset. 



initio molecular dynamics and phonon frequency calcula- 
tions. Ab-initio molecular dynamics calculations at 750 
K are performed with the same parameters as done for 
the monolayer ZnO in previous section indicate that the 
bilayer remains stable at 750 K after 3.0 ps. 



V. DIMENSIONALITY EFFECTS 

A comparison of Zn-0 bond length, cohesive energy, 
GGA bad gap, GWq corrected band gap calculated for 
3D wz-, zb-, h-ZnO, 2D bilayer and monolayer of ZnO 
are presented in Table-I, where interesting dimension- 
ality effects are deduced. These dimensionality effects 
are believed to be crucial for better understanding of 2D 
crystals. Quasi ID nanoribbons are not included to this 
discussion because of edge effects of the ribbon. 

Three dimensional crystals have larger number of near- 
est neighbors and also posses larger Madelung energy as 
compared to 2D crystals. wz-ZnO appears to correspond 
to the global minimum of ZnO II- VI compound. However 
the energy difference between wz-ZnO and 2D monolayer 
of ZnO is only ^0.5 eV and is smaller than one expects. 
It appears that sp 2 -like bonding which is stronger than 
sp 3 -like bonding and the 7r-bonding between adjacent p z 
orbitals, which contributes to stability by maintaining 
the planar geometry give rise to the relatively smaller 
energy difference between 3D and 2D monolayer ZnO. 
In the absence of significant 7r-bonding one expects that 
2D planar ZnO would be buckled (puckered) for stability 
as found in 2D Si and GaAs honeycomb structure! 12 * 13 ! 
Through puckered planar sp 2 orbital is dehybridized and 
is transformed to sp 3 -like hybrid orbitals. Here we note 
that population of oxygen sp 2 and p z orbitals is larger 
than that of corresponding Zn orbitals, since the former 
significantly higher electronegativity. 

Interestingly, while GGA bad gaps of 3D ZnO occurs 
0.75-0.96 eV, the band gaps of in 2D is 1.68-1.44 eV. The 
band gap in 2D is larger than 3D, since the energy differ- 
ence between sp 3 -like orbitals of Zn and O ions is smaller 
than the energy difference of sp 2 -orbitals. The latter 
gives rise to larger energy difference between bonding and 
antibonding orbitals. These band gaps are, however, un- 
derestimated by GGA, since they increase approximately 
three times upon GWq corrections. It appears that band 
gap values corrected by GWq is closer to experimental 
value and hence superior to GW correction. As for as 
dzn-o, Eq and Eq are concerned, 3D h-ZnO is inter- 
mediate between 3D and 2D structures. Also calculated 
values of bilayer ZnO are slightly closer to 3D than those 
of 2D monolayer of ZnO. 

Charge transfer from Zn to O is crucial for dimension- 
ality effects, but difficult to calculate precisely. Charge 
transfer, Sq, calculated using three different schemes, 
namely Lowdin, Bader and Siesta, yield different but 
consistent values. For example, the Lowdin values oc- 
cur around 1.4 electrons for 3D, but relatively smaller 
value of 1.35 electrons for 2D. Bader analysis yields 1.22- 
1.20 electrons for 3D, and relatively smaller value of 1.18 
electrons for 2D. As for Siesta, 3D values, they are 0.90 
electrons for 3D, but 0.87 electrons for 2D. Excluding 
the paradoxical situation with charge transfer values of 
zb-ZnO occurring close to those of 2D ZnO, the charge 
transfer between Zn and O is slightly smaller in 2D than 
3D. 
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ZnO 


Bonding 


dzn-O 


E c (eV/atom) 


GGA 


Go Wo 


GW 


GW 


Exp 


3D wz 


sp 3 


2.001-2.007 


8.934 


0.75 


2.76 


3.29 


2.4#3 


3.37^ 


3D zb 


sp 3 


2.001 


8.919 


0.65 


2.53 


3.09 


3.5EP 




3D h 


sp 2 +p z + VdW 


1.999 


8.802 


0.96 


2.84 


3.32 






2D bilayer 


sp 2 + p z + VdW 


1.932 


8.640 


1.44 


4.45 


5.10 






2D monolayer 


sp 2 + p z 


1.895 


8.478 


1.68 


4.87 


5.64 


3.5#31 





Table I: A comparison of bonding (or type of hybrid orbitals), the length of Zn-0 bond dz n -o, cohesive energy per Zn-0 pair 
Ec, band gap Eg calculated with GGA and corrected by Go Wo and GWo in the present study, earlier GW calculations and 
experimental value. 



VI. ZNO NANORIBBONS 

In this section, we consider bare and hydrogen pas- 
sivated armchair (a-ZnO) and zigzag (z-ZnO) ZnO 
nanoribbons. These nanoribbons are specified according 
to their width given in terms of n number of Zn-0 pairs 
in their unit cells. Hence, z-ZnO(n) indicates a zigzag 
ZnO nanoribbons having n Zn-0 pairs in their unit cell. 
We investigate their electronic, magnetic and mechanical 
properties. First, we start with the stability analysis of 
these nanoribbons. 



A. Stability Analysis 

We have analyzed the stability of the bare armchair 
ZnO nanoribbons having n = 9 Zn-0 atom pairs in the 
unit cell using the force constant method and the finite 
temperature molecular dynamics calculations. The dy- 
namical matrix was generated using the forces calculated 
in seven unit cells of a-ZnO(9). Results of this calculation 
were presented in Fig. [5] (a). Phonon dispersion curves of 
quasi ID nanoribbons, in general, show profiles expected 
from the folding of 2D phonon dispersion curves. Modes 
appearing above 700 cm -1 , however, are not expected 
from this folding. They were attributed to edge effects 
and reconstructions. 

Right panel of Fig. [5] (a) presents the acoustic region 
in ten times magnified scale. There are four acoustic 
modes, dispersion curves of which go to zero ask^O. To 
get exactly zero value, we have imposed the symmetries 
originating from translational and rotational invariance 
on the dynamical matrix. To impose these symmetries 
we have tuned the force constant matrix elements so that 
all forces on all atoms are zero when the whole structure 
is shifted in three dimensions or slightly rotated along 
the axis passing through the middle of the nanoribbon. 

Both the longitudinal and transverse acoustic modes 
have linear dispersions near the T point. Group veloc- 
ity of LA mode is slightly higher than that of TA mode. 
Out of plane ZA and twisting acoustical TW modes^l 
have quadratic dispersion near the T point, which is at- 
tributed to the rapid decay of the force constants with in- 
creasing neighbor distance. In fact there are some imag- 
inary frequencies in ZA mode near the T point. The ab- 



solute value of these frequencies do not exceed 0.5 cmT 1 
and they are purely an artifact of the precision of the 
numerical calculation. Using a finer mesh makes these 
imaginary frequencies disappear. 

We also have calculated vibrational modes of the fi- 
nite patch of ZnO having the length of five unit cells and 
the width of n = 9. Zigzag edges of this flake were sat- 
urated by hydrogen atoms to eliminate the magnetism 
and to simplify the numerical computations. Vibrational 
spectrum of this structure had no imaginary frequencies 
implying the fact that finite size flakes of 2D ZnO hon- 
eycomb structure is stable. The density of states, DOS 
of calculated phonon frequencies are presented in Fig. [5] 
(b). Note that, DOS calculated for three different honey- 
comb systems are similar. Owing to the edge effects the 
gap between the acoustical and optical branches of the 
ribbon is reduced. Similar effect as well as broadening of 
discrete mode frequencies cause DOS of the flake to de- 
viate significantly from that of 2D honeycomb structure 
in the gap. The vibrational modes attributed to strong 
Zn-H and O-H bonds are centered respectively at 1830 
cmT 1 and 3700 cmT 1 . These modes are not shown, since 
these frequencies are beyond the range of Fig. [5] (b). 

Ab initio molecular dynamic calculations were also car- 
ried out for a-ZnO(9) and a finite size patch of it. To 
eliminate constraints to be imposed by small unit cell, 
the infinite nanoribbon is treated by a supercell com- 
posed of five unit cells. Both nanoribbons are kept at 
750 K for 3.5 ps. As a result, infinite, periodic struc- 
ture of a-ZnO(9) composed of five unit cells had minor 
reconstructions at the edges, while its honeycomb struc- 
ture was preserved. Finite structure had the similar pat- 
tern at the armchair edges, but the reconstructions at the 
hydrogen saturated zigzag edges were more pronounced. 
These reconstructions made the whole system bend, but 
again the honeycomb structure was preserved around the 
central region of the nanoribbon. The results of this anal- 
ysis are interpreted that finite size ZnO nanoribbons are 
stable. 



B. Electronic and Magnetic Properties 

Bare and hydrogen terminated a-ZnO nanoribbons 
are nonmagnetic semiconductors with direct band gaps 
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Figure 6: (Color online) Atomic and electronic structure 
of bare and hydrogen terminated armchair nanoribbons a- 
ZnO(9). (a) Atomic structure with unit cell comprising n=9 
Zn-0 atom pairs, energy band structure and charge density 
isosurfaces of selected states of bare nanoribbons. The charge 
of the bands of edge states are localized at the edge atoms, (b) 
Same as (a) but for hydrogen passivated nanoribbon. Large 
medium and small balls indicate Zn, O and H atoms, respec- 
tively. Unit cells are delineated by dashed lines. 



which are relatively larger than that of the monolayer 
ZnO. Fig. [6] (a) and (b) shows the atomic and electronic 
structure of bare and hydrogen terminated a-ZnO(n) 
with n=9. The atoms at the edges of the bare a-ZnO are 
reconstructed; while one edge atom, Zn is lowering, adja- 
cent edge atom, O is raised. The energy band gap calcu- 
lated with GGA is direct and 1.92 eV wide and is larger 
than the band gap of 2D monolayer ZnO. Band decom- 
posed charge density analysis of a-ZnO(9) is presented 
in Fig. [6] The highest valance bands are degenerate and 
their charge accumulates on oxygen edge atoms for bare 
ribbon. The third band from the top of valence band 
is also associated with edge states. On the other hand, 
the charge of the lowest conduction band is distributed 
uniformly in the ribbon. The charge of the second con- 
duction band is mostly distributed at the dangling bonds 
on the edges of the ribbon. 

The passivation of Zn and O atoms at the edges by 
hydrogen atoms gives rise to changes in the electronic 
band structure. The bond lengths Zn-H and O-H bonds 
are calculated to be 1.53 and 0.97 A. As seen in Fig. [6] 
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Figure 7: (Color online) Variation of the calculated band gaps 
of bare and hydrogen passivated ZnO armchair nanoribbons 
with their widths n. Dashed line indicates the band gap of 
the infinite 2D ZnO. 



(b), upon passivation with hydrogen atoms, the recon- 
struction of edge atoms are weakened. At the end the 
edge state bands are discarded from the band gap and 
replaced by dispersive bands having charge distributed 
uniformly in the ribbon. The band gap of H-passivated 
a-ZnO(9) slightly increases to ~1.98 eV. 

The variation of the band gap Eq for the bare and hy- 
drogen saturated armchair nanoribbons as a function of n 
is given in Fig. [7] The band gaps are large for small n but 
approaches to those of 2D honeycomb structure as the 
width n — > co. This is an indication of the quantum size 
effect. For n<9, the value of the band gap of hydrogen 
passivated a-ZnO is significantly larger than that of bare 
ribbons; the difference practically disappears for n>20. 
In contrast to graphene and BN nanoribbons, family de- 
pendent variations of band gap^l is absent in Fig. [7] 
The variation of the band gap with n is an important 
property, which may lead to formation of quantum dot 
or multiple quantum wells through the size modulation P 

Earlier it has been 

reporteoPEH tnat a u a _z n O 
nanoribbons are semiconductors with a constant band 
gap of 2 eV. In our case, however, the band gaps of 
nanoribbons, especially for n<10, display apparent de- 
pendency on n. The tendency of the decrease in the band 
gaps to the 2D ZnO band gap is also observed when n is 
increased. 

In contrast to a-ZnO nanoribbons, all zigzag nanorib- 
bons (z-ZnO) we investigated in this work (namely 
n=5,7,9) show metallic character regardless of their 
widths. On the other hand, while all a-ZnO nanoribbons 
are nonmagnetic, bare z-ZnO nanoribbons are magnetic 
due to edge states. The magnetic properties of z-ZnO 
depend on whether the edge atoms are passivated with 
hydrogen. Our results regarding with the bare and hy- 
drogen passivated z-ZnO(n) nanoribbons for n=9 are pre- 
sented in Fig. [8] (a)-(c). In Fig. [8] (a), while the upper 
edge of the ribbon is made by O atoms, the lower edge 
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terminates with Zn atoms. The bare z-ZnO nanoribbon 
is magnetic, since the spin polarized calculations yield 
the total energy, which is energetically 34 meV favor- 
able than that obtained from spin unpolarized calcula- 
tions. The total magnetic moment of the system was 
calculated as 0.57 [Ib per unit cell. Figure [8] (b) shows 
total density of states (DOS) and band structure of bare 
z-ZnO(9) together with isosurface charge densities of se- 
lected bands. Bare z-ZnO(9) have isosurfaces of differ- 
ence charge density Ap occur around O edge atoms due 
to unpaired 0-2p orbitals. Clearly, bare z-ZnO is a fer- 
romagnetic metal. These results are in agreement with 
those predicting that the ferromagnetic behavior of ZnO 
nanoribbons due to un paired spins at the edges is domi- 
nated by oxygen atoms.E2I!H 

The band structure of the nanoribbon gives us further 
information about the magnetism of the system. When 
we plot the spin-up and spin-down bands we observe the 
splitting of these bands around Fermi level. The conduc- 
tion and lower valance bands are degenerate for spin-up 
and spin-down states, some of which are shown with yel- 
low isosurface charge densities in the figure. We also 
plotted charge density of two states of spin-split non- 
degenerate bands. Spin-down state is half-filled, while 
spin-up state is totally filled. Furthermore, their charges 
are confined to the O-edge of the ribbon. The spin polar- 
ization at the Fermi level is calculated in terms of den- 
sity of spin up and spin down states, I, E F ), namely 
P{E F ) = |D(t, E F ) -D(i, E F )\/[Dtf, E F ) + D(i,E F )} 
and is around 80 % for spin-down states although the 
spin-up states are majority spins. This suggests z- 
ZnO(n) with high spin polarization at E F can operate 
as spin valve device. 

As we discussed in Sec. |IV C| the magnetism of Zn- 
vacancy in the monolayer ZnO is directly related with 
O atoms having unpaired 2p orbitals. Similarly, bare 
z-ZnO nanoribbons have also magnetic states on O ter- 
minated edge of the ribbon. The splitting of spin-up 
and spin-down bands around Fermi level is removed when 
the edges of bare z-ZnO(9) in Fig. [8] (a) is passivated by 
hydrogen. The resulting structure is again metallic but 
nonmagnetic. The electronic properties of the z-ZnO(9) 
ribbon passivated with hydrogen is presented in Fig. [8] 
(c). Two bands just below the Fermi level is localized on 
the Zn edge of the ribbon and lowest conduction band is 
localized on the O edge of the ribbon. Similar results are 
also obtained for n=5 and n=7. The situation regard- 
ing the magnetism of bare z-ZnO is somewhat different 
from that of the bare zigzag graphene nanoribbons, which 
have ferromagnetic coupling along the edges, but antifer- 
romagnetic coupling between the edges P Moreover, un- 
like the case in z-ZnO nanoribbons, the magnetism of 
zigzag graphene nanoribbons are not destroyed upon ter- 
mination of edges with hydrogen atom. 

An important feature of of zigzag ZnO nanoribbons is 
that charged Zn and O atoms terminate different edges 
and thus induce sizable electric dipole moments. For 
the bare z-ZnO(9) the dipole moment is calculated to 



be 0.78 (electron x A) per cell. Upon H termination of Zn 
and O edge atoms the dipole moment increases to 1.07 
(electron x A) per cell and its direction is reversed. While 
the dipole effects are included in the electronic structure 
calculations of H-terminated z-ZnO(9) nanoribbons, the 
band gap under estimated by GGA may effect its metal- 
licity. Unfortunately, GWq corrections cannot be applied 
due to large number of atoms. 

C. Mechanical Properties of ZnO Nanoribbons 

The response of the nanoribbons to the applied uni- 
axial stress is crucial and provides information on the 
strength of the nanoribbons. More recently, the response 
of graphene to strain (or to tension) in the elastic defor- 
mation range has been an active field of study. For exam- 
ple, recent works^ indicated the effect of deformation on 
the electronic properties for band gap engineering. More- 
over theoretical studies^ have shown that carbon atomic 
chain can be derived from graphene under tension. C. Jin 
et al)^ showed that carbon atomic chains can be derived 
from graphene by electron irradiation inside a transmis- 
sion electron microscope. Therefore, the response of ZnO 
nanoribbons to tension is of crucial importance. 

Owing to ambiguities in defining the cross section of 
the ribbon one cannot determine the Young's modulus 
rigorously. We examined the variation of the strain en- 
ergy, E s (e) = E T (e) — E T (e = 0) with respect to the 
applied uniaxial strain, e = Ac/c, c being the lattice 
constant along the ribbon axis. The variation of E s (e) in- 
cludes information regarding the mechanical properties of 
the ribbon. For example, force constant, k ~ d 2 Er r jdc 2 
is obtained from the variation of E s . k is an observable 
and can specify the strengths of ribbons unambiguously. 
The effects of two edges due to incomplete hexagons be- 
comes important for narrow nanoribbons, but decrease 
with increasing width. In-plane stiffness C, which is in- 
dependent of the thickness, can also be used instead of 
Young's modulus. Defining Aq as the equilibrium area of 
the system, C can be given as 



By choosing a reasonable thickness "/i", Young's mod- 
ulus can be recalled as Y=C /h. The thickness value 
around 3.34 A was used in order to evaluate the elas- 
tic properties of SWNTs and graphene by some works 
in literature .H2HH] We note that while C is unique for 
2D infinite (periodic) honeycomb structure of ZnO, for 
nanoribbons it depends on how the width of the ribbon is 
taken in determining A . In fact, the width of the a-ZnO 
nanoribbon cannot be determined straightforwardly. It 
is not clear whether the distance is from Zn at one edge 
to Zn at the other edge, or from O at one edge to other 
O. One can also take the average of these two, or con- 
sider charge spill over from these edge atoms. Hence the 
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Figure 8: (Color online) Bare zigzag nanoribbon z-ZnO(n) with n=9: (a) Honeycomb structure with primitive unit cell 
delineated and difference charge density Ap of spin-up and spin-down states, (b) Total Density of spin-up and spin-down 
states (DOS), energy band structure and isosurface charge density of selected spin states. Large, medium and small filled balls 
indicate Zn, O and H atoms, respectively, (c) and (d) are same as (a) and (b). Yellow-light and blue-dark densities are for 
spin-up and spin-down states. Similarly red-continues and blue-dashed curves are spin-up and spin-down spin bands. Black 
lines are for nonmagnetic (spin-paired) bands. 



value of C is subject to small changes depending on how 
the width of the nanoribbon is taken. This uncertainty, 
however, diminish as the width increases; eventually the 
value of C converges to a fixed value. In this respect, 
once the width is fixed C, can be determined from k. 

We considered three nanoribbons, a-ZnO(9), a- 
ZnO(15) and z-ZnO(9). In order to eliminate the con- 
straints imposed by the periodic boundary conditions of 
small unit cell, calculations are performed using a su- 
percell composed of four primitive unit cells for a-ZnO 
and five unit cells for z-ZnO. The size of supercells is de- 
termined based on certain tests. Figure M (a) and (b) 
show the variation of the strain energy Ej[e) and band 
gap Eq{c) with the applied strain, e for both ribbons, 
a-ZnO(9) and a-ZnO(15). They display a parabolic E s 
vs e curves up to the strains e ~ 0.13. Beyond these 
strain values, energy vs strain variation becomes elastic 
but anharmonic. In the elastic range, the ribbons pre- 
serve their honeycomb like structure, but the hexagons 
are elongated uniformly along the direction of the strain. 
Zones of harmonic and anharmonic variation of E s (e) 



can be better deduced by calculating the tension force, 
F T (e) = -dE s (e)/dc. F T (e) is linear for < e <~ 0.13 
in the elastic-harmonic range, but becomes nonlinear for 
~ 0.13 < e as seen in Fig. [9]ja) and (b). If the applied 
tension is released, the elastic deformation in both har- 
monic and anharmonic ranges disappears and the ribbon 
returns to its original equilibrium geometry. 

The elastic deformation range of a-ZnO(9) and a- 
ZnO(15) ends with a sharp fall of the total energy. This 
point is identified as the yielding point of the ribbon oc- 
curring at e ~ 0.20. The yielding points are followed with 
structural transformation, where honeycomb cage struc- 
ture starts to undergo a change and part of the strain 
energy is relieved. Once the yielding point is passed the 
ribbon can stretch under significantly low tension and 
some fraction of the deformation will be permanent and 
irreversible if the tension is released. Similar behavior 
occurs when a nanowire of metals, such as Au or Cu; 
following an order-disorder transformation the elastic de- 
formation ends and the wire is elongated by one lattice 
constant (or it deforms plastically). However, in con- 
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Figure 9: (Color online) (a) Variation of the strain energy 
E a and band gap Eg of the bare a-ZnO(n) (n=9) with the 
applied uniaxial strain, e are shown by a curve with squares. 
Calculated tension shown by tin (red) curve is initially lin- 
ear, but deviates from linearity for large e. Atomic structure 
for three different values of strain. The boundaries of the su- 
percells comprising four unit cells and corresponding lattice 
constants c are indicated by vertical lines, (b) Same as (a) but 
for bare a-ZnO(n) (n=15). (c) Variation of force constant, ft 
and in-plane stiffness, C of a-ZnO with the width w of the 
nanoribbons. 



trast to the present case, the nanowire attains again the 
ordered state and start to deform elasticallyP 

The force constant is calculated to be k=2.05 and 3.88 
eV/A 2 for a-ZnO(9) and a-ZnO(15), respectively. The 
increase of k from 2.05 eV/A 2 to 3.88 eV/A 2 is due to 
the increase of the width from n=9 (~ 12.81 A) to n—15 
(~ 22.65 A). Note that if the Hooke's law were valid for 
these nanoribbons, the ratio «(15)//c(9) would be equal 
to the ratio of widths, namely 22.65/12.81=^1.77. The 
calculated value (3.88/2.05=1.89) is slightly higher. The 
discrepancy arises from the edges of the ribbons, which 
respond to tension differently from the central region. 
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Figure 10: (Color online) Variation of the strain energy E s , 
tension and magnetic moment /j, per supercell of the bare 
z-ZnO(n) (n=9) with the applied uniaxial strain. Atomic 
structures of the bare z-ZnO(9) for three different uniaxial 
strains. The supercells consist of five unit cells. 



Therefore, deviation from the Hooke's law becomes more 
serious as n decreases, but diminish as n — > oo. k val- 
ues of a-ZnO(9) and a-ZnO(15) calculated for primitive 
unit cell are 7.92 eV/A 2 (instead of 4x 2.05=8.20 due 
to numerical calculations performed in a large supercell) 
and 15.04 eV/A 2 (instead of 4x 3.88=15.52), respec- 
tively. These values are smaller than those corresponding 
to graphene ( 40.96 and 72.52 eV/A 2 ) and BN ( 34.02 
and 60.46 eV/A 2 ) nanoribbons. 

The calculated in-plane stiffness values C for a-ZnO(9) 
and a-ZnO(15) are 3.75, 3.71 eV/A 2 , respectivelyP The 
calculated C values of armchair graphene (BN) nanorib- 
bons for n=9 and n=15 are 18.45 (14.55) and 18.85 
(14.66) eV/A 2 , respectively. Graphene and BN nanorib- 
bons have significantly higher in-plane stiffness as com- 
pared to a-ZnO nanoribbons. The difference mainly orig- 
inate from the fact that a-ZnO honeycomb structure has 
larger unit cell than that of graphene and BN honeycomb 
structures. In view of the calculated n and C values, 
both graphene and BN nanoribbons are stiffer than ZnO 
nanoribbons. 

Figure [9] illustrates also the atomic structure of a- 
ZnO(9) and a-ZnO(15) nanoribbons for three different 
values of applied strain. After the yielding point, the 
honeycomb structure is destroyed and the polygons being 
smaller and larger than hexagons form. In some cases, a 
net magnetic moment can develop in those polygons. The 
onset of plastic range and variation of atomic structures 
can depend on the periodic boundary conditions. These 
artifacts of periodic boundary conditions are eliminated 
to some extent by using larger supercells comprising sev- 
eral unit cells. Nonetheless, the plastic deformations 
shown in Fig. [9] (a) and (b) have still large repeat period- 
icity. Even if the size and forms of polygons appear to be 
hysteric, their further investigation is of fundamental in- 
terest. The order-disorder structural transformation, as 
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well as dislocation generation appear to be absent in this 
study. The latter may require a treatment of deforma- 
tion by taking into account very large n and c. Moreover, 
as predicted for graphene nanoribbons, one expects that 
the onset of yielding in the presence of vacancy defects 
and also at high temperature can occur at small strain!^ 
Nevertheless, the plastic deformation in the present work 
is carried out under ideal conditions, hence its stochastic 
nature is not taken into account. The energy gap of the 
nanoribbon undergoes a significant change under uniax- 
ial tension. It decreases as e increases and becomes very 
small at the yielding point. Beyond the yielding point 
the band gap increases again. 

In Fig.[9](c) we plot the calculated values offeree con- 
stant, k relative to the width w of a-ZnO(n) taken as 
the distance between outermost atoms of both edges. As 
expected k is proportional to w and k(w) can be consid- 
ered as linear except small deviations from linearity due 
to edge effects. In fact, dn(w)/dw increases slightly with 
w. Whether a-ZnO(n) has even or odd n also causes to a 
small deviation from single linear variation. As expected, 
the variation of the in-plane stiffness C with w is not sig- 
nificant. However, C changes significantly between two 
subsequent values of n. 

The zigzag z-ZnO(9) nanoribbon displays a variation 
E s (e) similar to those of a-ZnO. The calculated value of 
k is 4.27 eV/A 2 , which is larger than that calculated for 
a-ZnO(9) having relatively smaller width. Similarly, in- 
plane stiffness value of z-ZnO(9) is slightly smaller than 
that of a-ZnO(9) and a-ZnO(15). C is calculated as 3.24 
eV/A 2 . The zigzag ZnO(9) remains to be ferromagnetic 
metal in the elastic deformation range, but their mag- 
netic moment decreases with increasing e. However, be- 
yond the yielding point it shows a sharp fall and diminish 
with increasing e in the plastic range as seen in Fig. [10] 
In contrast to a-ZnO, long polygons are aligned along the 
axis of the z-ZnO(9) in Fig. 
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VII. DISCUSSION AND CONCLUSIONS 

This work presents an extensive study on 2D mono- 
layer and bilayer ZnO, and ZnO nanoribbons together 
with their stability analysis. The monolayer of ZnO is 
an ionic and nonmagnetic, wide band gap semiconductor 
with significant charge transfer from zinc atoms to near- 
est oxygen atoms. Since DFT normally underestimates 
the calculated band gaps, the calculated band gap of spe- 
cific structures are corrected by GWq self-energy calcu- 
lations. ZnO has 2D hexagonal lattice forming a planar 
honeycomb structure. Our predictions, which contribute 
to a better understanding of this material are emphasized 
by way of conclusions: 

i) We have shown that, 2D ZnO monolayer and bilayer, 
quasi ID bare nanoribbons of ZnO and 0D patch of ZnO 
correspond to local minima on the Born-Openheimer sur- 
face and thus are predicted to be stable, ii) Ab-initio 
molecular dynamics calculations performed at high tem- 



perature corroborate our analysis obtained from the cal- 
culation of phonon frequencies, iii) We performed calcu- 
lations of phonon modes and density of states of phonon 
frequencies of 2D monolayer, quasi ID armchair nanorib- 
bon and a finite flake using force constant method. In 
particular, we revealed the acoustical twisting mode of 
armchair nanoribbon. iv) Our study of Zn-, 0-, Zn- 
O-vacancies and O-Zn antisite indicates that local mag- 
netic moments can be generated only by Zn-vacancies. 
v) We provided an extensive analysis of the electronic 
structure of hydrogen terminated, as well as bare arm- 
chair and zigzag nanoribbons. Armchair ZnO nanorib- 
bons are found to be nonmagnetic semiconductors. The 
band gaps vary with their widths. The narrow nanorib- 
bons have relatively larger band gap due to the quantum 
confinement effect. Bare zigzag ZnO nanoribbons are 
ferromagnetic metals. The atoms near the oxygen ter- 
minated edge of the ribbons acquire magnetic moments. 
The spin-polarization at the Fermi level may attain high 
values for specific zigzag nanoribbons. However, once 
O- and Zn- terminated edges are passivated with hydro- 
gen atoms, the zigzag nanoribbon becomes nonmagnetic 
metal. These electronic and magnetic properties might 
be useful for spintronic applications, vi) We found the 
minimum energy stacking of the ZnO bilayer, and pro- 
vided energetics and energy band structure correspond- 
ing to this stacking. Accordingly, bilayer ZnO is also a 
nonmagnetic, wide band gap semiconductor with slightly 
smaller band gap as compared to that of monolayer, vii) 
We provided an analysis of mechanical properties. ZnO 
nanoribbons under uniaxial strain show harmonic and 
anharmonic elastic deformation ranges and a yielding 
point. After yielding, the strain energy exhibits a sharp 
fall and the nanoribbon deform plastically. In the elastic 
range, hexagons are uniformly deformed, but honeycomb 
like atomic structure is maintained. After yielding point 
some of the hexagons are modified and reconstruct to 
different polygons which may be smaller or larger than 
hexagons. Variation of electronic and magnetic proper- 
ties with deformation and formation of polygons in the 
plastic deformation range are of fundamental interest. 
Calculation of force constants and in-plane stiffness in- 
dicate that the stiffness of ZnO nanoribbons is smaller 
than those of graphene and BN honeycomb structures. 

In summary, it is shown that single and bilayer ZnO 
and its nanoribbons in different orientations provide us 
for a variety of electronic and magnetic properties which 
may be interesting for further applications. Even if 
they have honeycomb structure common to monolayer 
graphene and BN, their properties exhibit important dif- 
ferences. 
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Through sequential order-disorder transformations a neck 
develops while the nanowire elongates. 
Because of uncertainties in measuring the width of a bare 
armchair nanoribbon, the value of C cannot be determined 
unambigiously. For example, by taking w as the distance 
between Zn atoms at both edges, w =12.186 A; alterna- 
tively iu=13.393 A if the width is taken as the distance be- 
tween O atoms at both edges of a-ZnO(9). The correspond- 
ing in-plane stiffness values are 3.977 eV/A 2 and 3.618 
eV/A 2 , respectively. One obtain C — 3.75 eV/A 2 . In con- 
clusions, while the value of C can be unambigiously deter- 
mined for 2D periodic ZnO honeycomb structure, the value 
of C calculated for a nanoribbon depends on its width. As 
w increases, C decreases, and eventually as w — > oo C 
approaches the value calculated for 2D periodic ZnO hon- 
eycomb structure. 



